Energy quantization and mean value inequalities 
for nonlinear boundary value problems 



We give a unified statement and proof of a class of wellknown mean 
value inequalities for nonnegative functions with a nonlinear bound on 
the Laplacian. We generalize these to domains with boundary, requir- 
ing a (possibly nonlinear) bound on the normal derivative at the bound- 
ary. These inequalities give rise to an energy quantization principle for 
sequences of solutions of boundary value problems that have bounded en- 
ergy and whose energy densities satisfy nonlinear bounds on the Laplacian 
and normal derivative: One obtains local uniform bounds on the comple- 
ment of finitely many points, where some minimum quantum of energy 
concentrates. 

1 Energy quantization 

One purpose of this note is to explain an 'energy quantization' principle that - in 
different forms - has successfully been applied to a variety of partial differential 
equations, such as minimal submanifolds, harmonic maps, pseudoholomorphic 
curves, and Yang-Mills connections. The common feature of these PDE's is an 
energy functional. (The solutions often but not necessarily are critical points 
thereof.) The phenomenon that we call 'energy quantization' is a consequence of 
a mean value inequality for the energy density. The second purpose of this note 
is to give a unified statement and proof of the underlying mean value inequality 
for the Laplace operator. In section 2 we moreover generalize this inequality to 
domains with boundary and inhomogeneous Neumann boundary conditions. 

In this section we generally consider a PDE for maps u : D — > T from a 
Riemannian manifold D (with possibly nonempty boundary 3D) to a target 
space T, e.g. another manifold, a Banach space, or a fibre bundle over D. The 
energy is given for all sufficiently regular maps u in the form 



where the integrand e(u) : D —> [0, oo) is a nonnegative energy density func- 
tion. Its key property is that for a solution u of the PDE, the positive definite 
Laplacian Ae(it) can be bounded above in terms of e(u) itself. 
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Abstract 




If this bound is linear in e(u), then one immediately obtains a C°-control on 
e(u) in terms of its mean values on geodesic balls. If moreover the energy E{ui) 
of a sequence of solutions Ui is bounded, then this provides a uniform bound on 
the energy densities e(u,) on any compact subdomain of D \ dD. In many cases 
this leads to a compactness result, i.e. to the convergence of a subsequence of 
the solutions Uj. 

For solutions of nonlinear PDE's however, the bound on Ae(zt) is usually 
nonlinear in e(u). In that case, for polynomial nonlinearities up to some order 
depending on the dimD, the mean value inequality only hold on geodesic balls 
with sufficiently small energy. In theorem 2.1 we give a precise and general 
statement of this wcllknown fact. In theorem 2.4, we then generalize this mean 
value inequality to domains D with boundary dD and bounds on the outer 
normal derivative -§^e\do- This allows to obtain uniform bounds on e{ui) up to 
the boundary for solutions Ui of a PDE with appropriate boundary conditions 
and with bounded energy. 

In the case of nonlinear bounds on the Laplacian or the normal derivative 
one only obtains locally uniform bounds on the complement of finitely many 
points: By a converse of the mean value inequality, the energy densities e(m) 
can only blow up at points where some nonzero quantum of energy concentrates. 
In the following lemma we give a blueprint for such energy quantization results. 

Here D is a Riemannian manifold (possibly noncompact or with boundary), 
A = d*d denotes the positive definite Laplace operator, and denotes the 
outer normal derivative at dD. For the sake of simplicity we make a technical 
assumption on the metric near the boundary. 

Assumption: A neighbourhood ofdD C D is locally isometric to the Euclidean 
half space H™. 

For a general metric near the boundary, the mean value inequality becomes 
harder and more technical to prove, but theorem 2.4 should generalize in the 
same way as theorem 2.1, and hence this lemma should extend to general Rie- 
mannian manifolds with boundary. 

Lemma 1.1 There exists a constant h > depending on n = dimD and some 
constants a,b>0 such that the following holds: Let a € C 2 (D, [0, oo)) be a se- 
quence of nonnegative functions such that for some constants A , A\, B , Bi > 

{ "+2 
Ae t < A + A x e, + ae t n , 
&\ gD ei <B + B iei + beT^. 

Moreover, suppose that there is a uniform bound I e { < E < oo. 

Jd 

Then there exist finitely many points, X\, . . . ,Xpj G D (with N < E/h) and 
a subsequence such that the are uniformly bounded on every compact subset 
of D \ {xi, . . . xjy}, and there is a concentration of energy h > at each Xj : 
For every 5 > there exists Ijj € N such that 

[ ei >h Vi> ij. s . (l) 



2 



Proof: Suppose that for some x € D there is no neighbourhood on which the ej 

are uniformly bounded. Then there exists a subsequence (again denoted (e,)) 

and D 3 Zi — > x such that ei(zi) = Rf with Ri — > oo. We can then apply the 

_ i_ 

mean value inequality theorem 2.1 on the balls Bg^Zi) of radius Si — R i 2 > 0. 
For sufficiently large i € N, these lie within appropriate coordinate charts of 
D. In case z e dD we use the Euclidean coordinate charts at the boundary to 
apply theorem 2.4, but we also denote the balls in half space by B$ i (zi). These 
mean value theorems provide uniform constants C and h := /x(a, b) > such 
that for every i G N either 

/ e 4 > h (2) 

JB 5 .( Zi ) 

or J Bs ei < ft, and hence 

R? = e(zi) < CA Sf + CB 6i + C(Af + B" + 5r n )h. 
In the latter case multiplication by 5™ = R i 2 implies 

Rf < caor;^ + cb r; ! ^ 1 + cn(AfR; q +b^r;^ + (3) 

As i — > 00, the left hand side diverges to oo, whereas the right hand side 
converges to Ch. Thus the alternative (2) must hold for all sufficiently large 
i S N. In particular, this implies the energy concentration (1) at Xj = x. 

Now we can go through the same argument for any other point x e D at 
which the present subsequence (ej) is not locally uniformly bounded. That way 
we iteratively find points Xj G D such that the energy concentration (1) holds 
for a further subsequence (a). Suppose this iteration yields ./V > E/h distinct 
points xi,. . . ,xn (and might not even terminate after that). Then we would 
have a subsequence (e^) for which at least energy h > concentrates near each 
xj. Since the points are distinct, this contradicts the energy bound J D ej < E. 
Hence this iteration must stop after at most \_E/h\ steps, when the present 
subsequence (e^) is locally uniformly bounded in the complement of the finitely 
many points, where we found the energy concentration before. □ 

The analogy in the compactness proofs for a variety of PDE's, including 
minimal submanifolds [A, CS] and harmonic maps of surfaces [SU], has already 
been observed and listed by Wolfson [Wo] . In all cases, the nonlinearities are ex- 
actly of the maximal order as in lemma 1.1. Here we discuss pseudoholomorphic 
curves and Yang-Mills connections in more detail. 

For pseudoholomorphic curves (with a 2-dimensional domain) the en- 
ergy is the L 2 -norm of the gradient, and the estimate Ae < C(e + e 2 ) leads to 
Gromov's compactness result [G, Wo]. For pseudoholomorphic curves with 
Lagrangian boundary conditions, these compactness results are also well- 
known. They can be proven via a specific choice of a metric for which J^e = 
(see [F] and [MS, 4.3]). Then the energy density can be extended across the 
boundary by reflection and the mean value inequality for 1" applies. For the 
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metric given by the almost complex structure, the Lagrangian boundary condi- 
tion only implies < C(e + e5 ), which fits nicely into our energy quantization 
principle. 

For Yang-Mills connections on 4-manifolds the energy is the L 2 -norm of 
the curvature. The bound Ae < C(e + e?) was used by Uhlenbeck [U] to prove 
a removable singularity result, which leads to Donaldson's compactification of 
the moduli space of anti-self-dual instantons [D]. For a proof of the energy 
quantization as in lemma 1.1 see also [We, Thm.2.1]. 

For anti-self-dual instantons with Lagrangian boundary conditions 
on H 2 x S, an argument along the lines of lemma 1.1 is used in [We, Thm.1.2]. 
There the energy density is given by the slicewise L 2 -norm of the curvature, 
e(A) = J s \Fa\ 2 - The Lagrangian boundary condition (which has global nature 
along the Riemann surface S) provides < C(e + e§), but one only has a 
linear bound Ae < ge with a function g that cannot be bounded in terms of e 
or a constant. The argument using (3) however allows for a mild blowup of A\, 
and the according estimate \g\ < CR 2 can be established. This result does not 
follow from the standard rescaling methods for Yang-Mills connections. 



2 Mean value inequalities 

In this section we state and prove the mean value inequalities that were desribed 
in section 1 and that the energy quantization principle is based on. 

The subsequent theorem is wellknown and proofs in an exhausting collection 
of cases can be found in the literature, e.g. [S, U]. Our aim here is to give a 
unified statement and proof. In theorem 2.4 below we moreover generalize this 
result to domains with boundary. 

We denote by B ro (x ) C R™ the open geodesic ball of radius r centered at 
x € R" and with respect to the present metric. The Laplace operator A = d*d 
and integration will also be using the metric given in the context. The Euclidean 
metric on R™ is denoted by its matrix 1. 

Theorem 2.1 For every neN there exist constants C , /i > 0, and 5 > such 
that the following holds for all metrics g on R" such that \\g — l||v^i,oo < 5. 

Let B r (0) C R™ be the geodesic ball of radius < r < 1. Suppose that 
e e C 2 (£> r (0), [0, oo)) satisfies for some Aq, A\, a > 

n+2 f _„ 

Ae<A + Aie + ae « and / e < [ia 2 . 

Js r (o) 

Then 

e(0) < CA Q r 2 + C(AJ + r~ n ) [ e. 

JB r {Q) 

Remark 2.2 By using local geodesic coordinates the above theorem also implies 
a mean value inequality on closed Riemannian manifolds with uniform constants 
C, [i > 0, and for all geodesic balls of radius less than a uniform constant. 
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The special case Aq = A\ = a = of theorem 2.1 and the starting point for 
the proof is Morrcy's [M] mean value inequality for subharmonic functions. A 
proof of the version below can be found in e.g. [LS, Thm.2.1]. For the Euclidean 
metric g = 1 we give an elementary proof in lemma 2.5 below. 

Proposition 2.3 For every n £ N there exist constants C and S > such 
that the following holds for all < r < 1 and all metrics g on R n with 
\\g - l||iyi,oo < 6: Ife£ C 2 (B r (0), [0,oo)) satisfies Ae < then 

e(0) < Cr- n [ e. 

Proof of theorem 2.1: 

This proof is based on the Heinz trick, which is to consider the maximum c 
of the function / below. This allows to replace the bound on the Laplacian 
by a constant depending on c. One then obtains the result from the mean 
value inequality for subharmonic functions and a number of rearrangements in 
different cases. 

Consider the function f(p) = (1 — p) n sup Bpr ( ) e for p £ [0, 1]. It attains its 
maximum at some p < 1. Let c = sup B _ r ( ) e = e(x) and e = ^(1 — p) < |, then 

e(0) - /(0) < f(p) = 2"e"c. 

Moreover, we have for all x £ B sr (x) C B r {y) 

e(x) < sup e = (l-p-e)- n f(p + e) < 2"(1 - = 2"c, 

B (p+£)r(0) 

and hence Ae < A a + 2 n Aic + 2 n+2 ac i ^ 2 . Now define the function 

v(x) := e(x) + ±(A + 2 n c(A 1 + 4aci))\x - x\ 2 

with the Euclidean norm \x — x\. It is nonnegative and subharmonic on B er {x) 
if the metric is sufficiently C 1 -close to 1. This is since Aj|x — x\ 2 = 2n for the 
Euclidean metric and |x — x\ < er < 1 is bounded, so A\x — x\ > n whenever 
|,g — l|| w i,oo < 5 is sufficiently small. The control of the metric also ensures 
that the integral J B ^ \x — x\ 2 is bounded by the following integral over the 

Euclidean ball B% r {x) 

2 f \x-x\ 2 = 2 f P t" +1 Vol5"- 1 dt = C n (pr) n+2 . 
JbI pt (x) Jo 

So we obtain from proposition 2.3 with a uniform constant C for all < p < e 

c = v(x) < C{pr)- n { v < C(A +c(A 1 +4aci))(pr) 2 + C(pr)- n j e. 

JB pT (x) J B pr {x) 

If CA (sr) 2 < |c, then we can drop A from this inequality, just changing the 
constant C. Otherwise e(0) < c < CA r 2 already proves the assertion. 
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Next, if C(A 1 + Aaci)(pr) 2 < \, then this implies c < 2C(pr)~ n J Br{Q) e. So 
if C(A\ + 4ac")(er) 2 < | then p = e proves the assertion, 

e(0) < 2" £ "c < 2 n+1 Cr- n J BA0) e. 

2 — - 

Otherwise we can choose < p < e such that pr = (2C(A\ + 4ac~)) 2 to 
obtain with a new uniform constant C 

e(0) < c < C(A 1 +Aaci)^ J BpAs) e. 

Again we have to distinguish two cases: Firstly, if 4ac« < A\ then this yields 

e(0) < C(2A^ j BpAs) e. 

Secondly, if A\ < Aaci then c < cC(8a)% J B ^ e and thus with some p > 

Js r (o) e > ^ a_t - 

So altogether we either have the above or with some uniform constant C 
e(0) < CA a r' + C(Af+r- n )J Br(0) e. 

□ 

Next we generalize the mean value inequality theorem 2.1 to the half space 

H" = {(x ,x) | x e [0,oo),z e R"" 1 }- 

In order to generalize the mean value inequality to manifolds with boundary 
we would have to consider general metrics on H™. This would however disturb 
the elementary geometric proof, so we restrict this exposition to the case of 
the Euclidean metric. In that case the outer normal derivative ^\du n is just 
— g§^Uo=o- We moreover denote the intersection of a Euclidean ball with the 
half space by 

D ro (x ) := B ro (x ) ("1 H". 

Theorem 2.4 For every n > 2 there exists a constant C and for all a, b > 
there exists p(a, b) > such that the following holds: 

Let D r (y) C H™ &e the Euclidean n-ball of radius r > and center y € H™ 
intersected with the half space. Suppose that e e C 2 (D r (y), [0, oo)) satisfies for 
some A , Ai, B , B x > 

/ Ae < A) + Aie + ae^, / . . 

< a, n+i and / e<p(a,b). 

[ 9u| aH „e <B + B ie + be - , 

TTien 

e(y) < CA r 2 + C5 r + C(Af + + r~ n ) f e. 

JD r (y) 
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We will prove this theorem in three steps. The first step is the generalization 
of proposition 2.3 to domains with boundary and subharmonic functions in the 
sense of the weak Neumann equation: A distribution e on a manifold M is called 
subharmonic if for all ip e C°°(M, [0,oo)) with |^| 9M = 



> / eAtp = / -0Ae + 

' M J M 



JdM 



For e G C 2 (M) the equality above holds and implies that Ae < and -^\qm < 0. 

Lemma 2.5 There exists a constant C such that the following holds for all 
R > and y G H": Suppose that e G C 2 (D r (y), [0, co)) satisfies 

Ae < 0, 
Then 

e(y) < Cr- n [ e. 

JD r (y) 

Proof: We write H" = {(x a ,x) | x G [0, oo),x G R n_1 } and use spherical 
coordinates (xo, x) = (yo + r cos0 , y + r sin • z) =: (r, 0, z) with r G [0, oo), 
G [0,7r], and z G S" 1 " 2 C K". (For n = 2 this notation means 5° = {-1, 1}, 
and integration j s0 . . . dvol s o will denote summation of the values at these two 
points.) Now the boundary of D r (y) has two parts, 

Z r := dD r (y)ndW l = {(0 , x) \ |x - y\ 2 < r 2 - y 2 }, 
r r := dB r (y) DM n = {(r,</>,z) I G [0,Mr)],z G S n - 2 }. 

Here we use <j>o(r) : = arccos(— yo/r). For y n > r we set ^oM := 7r, so T r is the 
entire sphere and the set Z r is empty. With this we calculate for all r > 



d_ 

r 

r4>o(r) 



00 (f) 



i . r ™ +1 / / e(r, 0, z) (rsin^)™ 2 dvol s „-2r< 



r<t>o(r) r 

= / 9 r e(r, 0, z) (sin0)" -2 dvol s „-2 d</> 

Jo is"- 2 

+ ^ ^ 2e (r,0 o (r),z)(sin0o(r)r- 2 dvol s „- 2 . 



(4) 



Note that (f>o(r) is constant for y — as well as for r < y . So firstly in case 
yo > we have for all < r < y 



d 
dr 



[ e] = r- ,l+1 f £e = - r - ,l+1 f Ae > 0. (5) 

V JT T J JdD r {y) " JD r {y) 
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In that case we moreover have 



lim (r~ n+1 f ej = Vol S"^ 1 e(2/), 
so integrating J Q T r n_1 . . . dr proves the lemma for all R < 2y , 

R 

L2- n R n \o\S n - 1 e{y) < f* f edr < [ e. 

Jo Jr r J Dr (y) 



(6) 



Next, in case yo = we have for all r > 
d 



dr 



I Ae - ,-»+' / &e > 0. 

JD r (y) JZ r 



= _ r ~n+l 



Since lim r ^ (r n+1 J r ej = iVolS 1 " 1 e(y), integration over < r < R then 
proves the lemma for y = and all R > 0, 



^iTVolS"" 1 ^) < I" f edr = [ 

JO JT r JL 



e. 

D R (y) 



Finally, in case R > 2y n > we obtain from (4) for all r > y 

d f.-n+l f c \ ^ ~V0 f , , ^(■j,,^n-2 



r-+M e > / e(r,Mr),z)(sm<P (r)) n - 2 dvol s ^. 

dr \ Jr r J ryjr 2 - y 2 Js™- 2 

Now we can use (6), (5), and integrate the above to obtain for all yo < r < \R 
VolS n - 1 e(y) <r- n+1 f e 

+ [ yoP^ n (p 2 -Vo)^ I e(p,Mp),z) dvol s „-2 dp. 

Jyo JS n - 2 

Since (p, (f>o(p), z) € dW 1 , we already know that 

2n f 2 n+1 n f 

^MPU) < e < e . 

2 

With this we find that for all < y < r < \R 



S 71 - 1 e(y) < r- n+1 [ e 

+ - - ; n — 7- / t- 2 (l-t- 2 ) 2 dt \ 
< r- n+1 f e + C n R- n f 



'D R (y) 



Dr(v) 

(7) 



Here we have introduced a constant C n that only depends on n > 2, in particular 
on the value of the integral in t: For n = 2 we calculate it explicitly, 

/rVo 1 _i r i ry^ 1 

r 2 (l-r 2 ) 2 dt = |arccos(t _1 )J = arccos(^) < f. 



For n > 3 we have 

n-3 /""fo 



f r V n-3 f r Vo 

J r 2 (l-t- 2 )^ dt < j i^dt = 1 - £ < 1. 

Now from (5) we know that (7) also holds for < r < yo (with C n = 0), so 

r — 

integrating J Q 2 r n ~ Y . . . dr proves the lemma in this last case, 



IdfVol^eO/) < I 2 I edr + I(f) n C n ir" f e < C f 
Jo Jr r Jdr(v) Jd 



e. 

io Jv r Jd r (v) Jd r (v) 

□ 

Proof of theorem 2.4: 

With lemma 2.5 in hand, the second step of the proof is to assume constant 
positive bounds, Ae < A and J^| aH „e < B and find a constant C such that for 
all r > and y e H™ 

e(y) < Cr- n [ e + CAr 2 + CBr. (8) 
JdAv) 

That is, we first prove the theorem with A\ = B\ = a = b = 0. To do this 
consider the function 

v(x) := e(x) + ±A\x - y\ 2 + (B + ±Ay )x . 
It is positive and satisfies Av < and -§^\ m n v < 0, so lemma 2.5 implies that 

e(y) = v(y) - (B + ±Ay )y < v(y) < Cr~ n f v. (9) 

JD r (y) 

In case r < y we just use v(x) = e(x) + j^A\x — y\ 2 , then the same holds, and 
moreover 



/ - = I e + ft^YolS^dt = j e+ ^g^4f 

JD r (y) J D r (y) JO J ' D r (y) 

In case r > yo we have (using Xo < 2r on B r (y)) 

[ v < [ e + ±A f r+'VolS^dt + (B+±Ay ) f 

JD r (y) J D r (y) JO J D r 

< ( e + v °\ S Zi\ Ar n+2 + (B + ±Ar)?-Vo\S n - 1 r n+1 

— / 2n(n+2) v n ' n 

JD r {y) 

In any case, putting this into (9) proves (8). 



n+2 
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Finally, to prove the theorem we consider - analogous to the proof of theo- 
rem 2.1 - the function f(p) = (1 — p) n su\) DpAv ^ e defined for p e [0, 1]. It attains 
its maximum at some p < 1. We denote c = sup^^) e = e(x) and e = \{l — p), 
then e(y) < 2 n s n c and e(x) < 2 n c for all x e £> £r "(z). Thus on D er (x) C £> r (y) 
we have Ae < A + 2 n c{A 1 + Aaci) and £\ mn e < B 1 + 2 n c(B 1 + 2bc^). 
Putting this into (8) yields for all < p < e 

c = e(x) < C(pr)- n J DpAs) e + C(A + 2 n c(A 1 + 4ac*)) (prf (10) 

+ C(B + 2 n c(B 1 +2bci))pr. 

To deduce the claimed mean value inequality from this, we have to go through a 
number of different cases. Firstly, if CA (er) 2 + CB er > ±c, then since e < i 

e{y) < c < CA r 2 +CB r, 

which proves the theorem. Otherwise (10) continues to hold with Aq and Bq 
dropped (and another constant). Next, let < e' < e be the solution of 
A^e'rf + Bxe'r = 2- ,l - 1 C'- 1 or in case A^erf+B^r < 2-™- 1 C"- 1 let e = s' . 
Then we can rearrange (10) to obtain for all < p < e' and yet another constant 

J DAy) e > J DpAs) e > c(prr(C-i - aci{prf bc^pr). (11) 

Now if a, b ^ let 77 (a, b) > be the solution of 

ar 1 2 + br ] = \C~ X . 

If C" pr = rj(a, b) for some < p < e', then the theorem holds by the following 
definition of p(a, b) > 0, 

lD r (y) e ^ ib^M)" = : MM)- 

Otherwise we must have cMe'r < T](a, b), so (11) with p — e' gives 

c < 2C(e'r)- n J Dr{y) e. (12) 

In the special case a = b = we get the same directly from (11). In case e' = e 
this proves the theorem since e(y) < 2 n e n c. Otherwise e' < e satisfies with 
another constant 

2C- 1 = A^e'r) 2 + B x e'r + CT 1 

= {^/A~ 1 e l r + C-^) 2 + {B 1 -2C-^^/A' 1 )e'r 

= (^B l£ 'r + C-i) 2 + (A 1 - %B 2 )(e'r) 2 . 

From this one sees that either B 1 < 2C~^^/A^ and e'r > (\/2 - l)C~iA 1 2 
from the second line, or A\ < ^B\ and e'r > 2{\f2 — l)C~ 1 B^ 1 from the third 
line. Putting this into (12) we finally obtain in this last case 

e(y) < c < C(Af + B?) J DAy) e. n 
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